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We use Lorentz invariance and the QCD equations of motion to study the evolution of 
functions that appear at leading order in a 1/Q expansion in azimuthal asymmetries. This 
' includes the evolution equation of the Collins fragmentation function. The moments of these 

, functions are matrix elements of known twist two and twist three operators. We present the 

^SJ ' evolution in the large Nc limit, restricting to non-singlet for the chiral-even functions. 
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CN . I. INTRODUCTION 

Azimuthal asymmetries in hard scattering processes with at least two relevant hadrons constitute a rich 
phenomenology, suitable for studying quark and gluon correlations in hadrons. By relevant hadrons we mean 
t — ■ hadrons used as target or detected in the final state. A well-known azimuthal asymmetry appears in the semi- 
' inclusive deep inelastic polarized leptoproduction of pions {ep^ — > e'lrX) generated by the so-called Collins 
. effect [0 . This asymmetry is one of the few possibilities to gain access to the so-called transversity or transverse 
spin distribution function j^,^ , which is the third distribution function needed for the complete characterization 
of the (coUinear) spin state of a proton as probed in hard scattering processes. In contrast to the transversity 
function, the evolution of the Collins fragmentation function has not been presented yet. Knowledge of this 
[ evolution is indispensable for relating measurements at different energies. 

In processes like the semi-inclusive leptoproduction mentioned above, it is important to take transverse 
O i' momentum of partons into account. The parton distribution functions as a function of a light-cone momentum 
D , fraction and transverse momentum have first been studied by Ralston and Soper ||^ for the Drell-Yan process at 
' tree level. The precise operator definition of such transverse momentum (p^) dependent distribution (and also 
fragmentation) functions is a non-trivial issue (mainly because of gauge invariance) and in several studies this 
matter has been addressed [Q-Q. Besides their definition and appearance in cross sections, they have not yet 
. , been considered beyond tree level. One thing that one wants to know is how these Pt dependent distribution 
' and fragmentation functions evolve, for instance the fragmentation function pointed out by Collins, which is 
one of the new functions entering in the description of hadrons when transverse momentum is considered. Its 
evolution will be one of the new results presented in this paper, although we limit ourselves to the large Nc 
limit, in which case the evolution is an autonomous one. 

In general, factorization of hard scattering processes means that it is possible to separate parts containing 
only soft or hard physics. In the perturbative calculation of the hard subprocess (the partonic cross section) one 
encounters collinear divergences, which one can absorb into the matrix elements (or equivalently the distribution 
or fragmentation functions) describing the soft part of the process. This redefinition of the soft parts should be 
possible to all orders in the coupling constant. This procedure introduces a factorization scale and the goal is 
to calculate the dependence on this scale, which determines the high-energy behavior of the cross section. 

Factorization crucially depends on the presence of a large energy scale in the process, such as the space-like 
momentum transfer squared in leptoproduction or the time-like momentum squared of a lepton pair in Drell- 
Yan scattering. In this paper we will be concerned with functions that appear in processes which have, apart 
from such a hard scale, an additional soft momentum scale, related to the transverse momentum of the partons. 
The first factorization theorem for such a situation was obtained in Ref. [|j for the process e~^e~ hih2X, 
where the vector boson has a large invariant mass Q, but a small transverse momentum Qj, with respect to the 
two almost back-to-back hadrons hi and h2, i.e. ^ Q^. Similar situations occur in the Drell-Yan process, 
where one has besides the momentum of the lepton pair, two hadron momenta and in one-hadron inclusive 
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leptoproduction where one also deals with three momenta: the large momentum transfer, the target momentum 
and the momentum of the produced hadron. 

The effects of parton transverse momenta lead to azimuthal asymmetries in such processes, often coupled to 
the spin of the partons and/or hadrons. Just as for spin asymmetries, the azimuthal asymmetries provide a 
rich new phenomenology in Drell-Yan scattering, semi-inclusive deep inelastic scattering and electron-positron 
annihilation |p|,^-p^ . In this paper we will study the scale dependence of the various distribution and frag- 
mentation functions appearing in these (polarized) processes. We do this for specific moments in both Pt and 
X, employing Lorentz invariance and the QCD equations of motion. The moments in x for leading (coUinear) 
distribution functions (appearing for instance in inclusive leptoproduction) are related to matrix elements of 
twist two operators. On the other hand, for the transverse moments entering the azimuthal asymmetry expres- 
sions of interest, one finds relations to matrix elements of twist two and twist three operators, for which the 
evolution, however, is known. In the large Nc limit this evolution becomes particularly simple and is known 
from studies of inclusive processes where the transverse moments can be eliminated from the expressions. This 
knowledge allows us to obtain the evolution equations for the desired Pt and x moments of the transverse mo- 
mentum dependent distribution and fragmentation functions, that enter in azimuthal (spin) asymmetries. Such 
asymmetries recently have gained in interest, as can be seen from the experimental studies in Refs. |Ij-Eq|. 



II. FORMALISM 



In this paper we will study the scale dependence of the distribution and fragmentation functions in (for 
example) one-hadron inclusive leptoproduction {eH — > ehX) at leading order in an expansion in where = 
—Q^ is the space- like momentum transfer squared. Experimentally, we are interested in azimuthal asymmetries 
in the current fragmentation region, in which case the target hadron momentum P and the produced hadron 
momentum Ph satisfy P-Pt ~ Q^- We introduce two light-like vectors via the hadron momenta and parameterize 

,= -^n_--^n+ + ,., (3) 

where n+ and n_ are two light-like vectors, chosen such that n+ • n_ — 1. We will often refer to the ± 
components of a vector p, which are defined as p^ = p ■ n^. We define the transverse momentum scale 
= — gj. We are interested in the region where Q\. ^ . Up to mass corrections of order one then 

has = — Q"^ ~ Q^. The ratio ^ — —q'^/P^ w Q^/2P ■ q — Xb is the Bjorken scahng variable and the 
ratio C ~ Ph ^ P ■ Ph / P ■ q — Zfi is the usual fragmentation variable. 

In the case of inclusive deep inelastic scattering the soft part of the process is described by a correlation 
function. To be more specific, in leading order in powers of 1/Q, the forward scattering Compton amplitude 
T^'^ can be written as 

Tf"" = J dxTi [S'"'{x)^x)] + ■■■ (4) 

Here, S*^" and <I> are the hard and soft scattering parts, respectively. The hard part S^'^ is the 7*-quark forward 
scattering amplitude. The color gauge invariant soft part is defined as 



$,,(x)^ I ^ e^P<{P,S\^^{0)U{0,OMO\P,S) 



(5) 

LC 



where the subscript 'LC indicates — — and 



U{0, 0=V exp l^-ig ^ dr^- A+ (r;) j , (6) 

is a gauge link with the path running along n_. The correlation function $ is a function of the light-cone 
momentum fraction x = p^ /P^ only. The most general parameterization of $ which is in accordance with the 
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required symmetries (hermiticity, parity, time reversal) and which is relevant for calculations at leading order 
in 1/Q is given by 

$twi.t~2(^) ^ 1 {/i(a;)^+ + ^,gi(x) 75 7i+ + /ii(a;)75^.^+} , (7) 

where also for the spin vector a decomposition in n± is adopted, S — {P^ /M) n+ — {P^ /M) n_ + St- 
Specifying also the flavor one also encounters the notations q{x) = fi{x), Aq(x) = gfix) and Sq{x) = ATq{x) = 
hl{x). The evolution equations for these functions are known to next-to-leading order and for the singlet 

/i and gi there is mixing with the unpolarized and polarized gluon distribution functions g{x) and Ag{x)^ 
respectively. 

Denoting these functions as twist-2 makes sense because the local operators connected to the McUin moments 
of these functions are related to the matrix elements of local twist-2 operators, like '0 7^ (£>+)" ip. 

If one calculates T^^ up to order 1/Q one needs also the M/P^ parts in the parameterization of ^{x) 



^tw.st-.(^) = ^ e{x) + grix) 75^. + h^x) 75 

iSi^ eL(a:)75 - frix) e^^jpSTa +ih(x) \ . (8) 



2P- 

We have not imposed time-reversal invariance in order to study also the T-odd functions, which are important 
e.g. in the study of fragmentation functions. The fmictions e, gx and are T-even, the functions e^, Jt and 
h are T-odd (we will not concern ourselves with the formal problems related to T-odd distribution functions 
[0,0). The leading order evolution of e, gr and is known |l9|-p^ and for the non-singlet case this also 
provides the evolution of the T-odd functions cl , /t and h respectively, for which the operators involved differ 
only from those of the T-even functions by a 75 matrix. 

The twist assignments of these functions is better seen by considering the light-cone correlators 



$S,,(x,y)^ ^ / ^^e^^-(«-'')+*^ ''(P,^|^,(0)Z^(0,r;)*7^?(,7)W('?,O^.(OI^,5) 



(9) 
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depending on a; = / P^ and y = / P^ . It is parameterized in terms of two-argument functions 
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+ Hd (x, y) S, 757? ^+ + Ed{x, y) 7? 1 , (10) 

where parity invariance has been imposed. Hermiticity leads to G*j^{x,y) = —GDiy,x), G^{x,y) — GD{y,x), 
H^(x, y) = Huly, x), and E'^{x, y) — —Eoiy, x). Time reversal invariance would require these functions to be 
real. The QCD equations of motion can be used to relate the twist-3 functions appearing in the parameterization 
of $ to the correlators <I>^. This gives the relations 

dy Goix^y) + GD{x,y) = x gxix) - —hi{x) + ix frix), (11) 

/Tfl 
dy HD{x,y) =xhL{x) - gi{x) ~ ixeL{x), (12) 

dy ED{x,y) ^xe{x) ~ — fi{x) + ixh{x). (13) 

The local operator matrix elements corresponding with the moments of the functions in $^ (x) are (note that 
a is transverse) twist-3 operators, up to quark mass contributions multiplying twist-2 operators. 



III. TRANSVERSE MOMENTUM DEPENDENT DISTRIBUTION FUNCTIONS 

If one considers a semi-inclusive hard scattering process in which two hadrons are identified (in either initial 
and/or final state), then the treatment of transverse momentum is important. One needs to study correlation 
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functions that also depend on the transverse momentum, ^{x,Pt), for which the most general parameterization 
involves more functions. To be more specific, one needs the lightfront correlation function 



(27r)3 



e^P<{P, S\ ^MU{Q, -cx3)Z^(-oo,0 ^.(0|P, S) 



(14) 



At leading orders in powers of 1/Q the foUowingtransverse momentum dependent distribution functions are 
needed to parameterize this correlation function [p|,pl]| 



M 



5is(a;,Pr) ji+ls - hiT{x,pl) ia^^j^Sl^n'^^ 



htix,PT) 



M 



+ ht{x,pi.) 



2 \ ^ ^yPT^-\- 



M 



(15) 



We used the shorthand notation 

gis{x,pT) = Sl giL{x,p: 



2 N I [Pt ■ ) / „2 N 

^ ' — — giT{x,p^), 



(16) 



and similarly for h^^. The parameterization contains two T-odd functions, the Sivcrs function |p3| , p4| and 
the function hj, the distribution function analogue of the Collins fragmentation function ^ . Upon integration 
over Pt, Eq. (|l5|) reduces to <i>(x) with 

(17) 
(18) 
(19) 



/i(x) = J (fpj,fi{x,pl), 
gi{x) = J (fpT,giL{x,pl) 
hi{x) = / (fpT 



hiT{x,pl) + hir^{x,pl) 

The deep inelastic scattering process is only sensitive to the latter three functions, but in semi-inclusive deep 
inelastic scattering or in the Drell-Yan process (at small g^), one is sensitive to the p^-dependent functions. 
At measured Qj- one deals with a convolution of two p^'dependent functions, where the transverse momenta of 
the partons from different hadrons combine to Qt P, [io| , ^5|] . A decoupling is achieved by studying cross sections 
weighted with the momentum g", leaving only the directional (azimuthal) dependence. The functions that 
appear in that case are contained in 



which projects out the functions in (f>(x,pr) where appears linearly. 



(20) 



-g^^ix)S^Tri,,,-S.ht^^\x)^:^^^ 



ft^'\-)e\.,rn'^S^-^ht^^\x)^A^ 



(21) 



where we define p^/2M2-moments (transverse moments) as 



f pI 



V2M2 



f{x,PT)- 



(22) 



The functions hj- and are T-odd, vanishing if T-reversal invariance can be applied to the matrix element. 
At this point one can invoke Lorentz invariance as a possibility to rewrite some functions. All functions in 
and ^gix) involve nonlocal matrix elements of two quark fields. Before constraining the matrix elements 
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to the light-cone or lightfront only a limited number of amplitudes can be written down |10) . This leads to the 
following Lorentz-invariance relations il9|,|io|,0 

(23) 
(24) 
(25) 
(26) 

From these relations, it is clear that the p^/2A/^ moments of the p^-dependent functions, appearing in $S;(x), 
involve both twist-2 and twist-3 operators. 

Starting from the defining expression of one obtains, after weighting with p^, the gauge invariant operator 
structure, 



yi 




flL = 




fT = 


d ±(1) 


h = - 





i+=iT=o 



(27) 



To see the partonic interpretation of the second term, consider the light-cone gauge = 0), in which case 
= d-A" and the gauge links become unity. Up to some for our purposes not relevant complications with 
boundary terms, the expression thus reduces to 



mx) = J dy [^"o{x,y)-n{^,y)] 



(28) 



where the leading part of the <i>^ matrix element (in gauge = 0) is built from ipj^ and A" fields. The correlator 
<i>^ can be parameterized analogous to $^ with (interaction-dependent) functions Ga, Ga, Ha and Ea with 
similar hermiticity properties as the functions in i^j- Using Eq. (^) we define the following combinations. 



dy 



GA{x,y) + GA{x,y) ^ x grix) ~ ^ hi{x) - g[l^{x) + i x fT{x) + f^rj,^\x) = x grix) + ix frix), (29) 



2 I dy HA{x,y) ^ xKl^x) - ^ gi{x) +2h'^l^\x) - ixeL{x) = xhL{x) - ixeL{x), 
2 j dy EA{x,y) = X e{x) — ^ fi{x) + i x h{x) + 2 hj^^^\x) = x e{x) + ix h{x) . 



(30) 
(31) 



In principle, one can connect the functions defined here to those appearing in the treatments of Ellis, Furmanski 
and Petronzio j2^ or to those in the treatment of Jaffe and Soldate [^. We end this section by giving the 
relation to the functions used in a more recent treatment by Belitsk y [p8| , p9| ; comparison of the equations of 
motion Eqs. ( p9| ) - ( pT| ) and Lorentz invariance relations Eqs. ( |2^ ) - (|26|) with those given in Ref. |^ leads us 
to identify 



g[}^ix)=Kix) 



dyfiy), 



with K, D as defined in Ref. 



xgrix) = J dx'D{x,x') + J dyf{y), 
and 

,D{y,x') + Dix',y) 



fiy) 



dx' 



x' -y 



(32) 
(33) 

(34) 
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IV. RELATIONS BETWEEN TWIST-3 FUNCTIONS AND TRANSVERSE MOMENTS 



Using the equations of motion relations in Eqs. ( |29| ) - (^ij) and the relations based on Lorentz invariance 
in Eqs. ( p3| ) - (p6|), it is straightforward to relate the various twist-3 functions and the p'^/2AP (transverse) 
moments of p^-dependent functions. The results, grouping relevant combinations, are 



9t{x) = 
sir (2;) 



y 

9i{v) 



m 
M 



hi(x) 



dy 



y 



m 
M 



X 

\ hi{y) 
dy — — 



\ hi{y) 
dy — ^ 



dy 



9t{x) 

griy) 
y 



(35) 
(36) 



hL(x) = 2x I dy 1 



M 



^-2x f\y''^y^ 



ht^^\x) 



hi{y) m f gi{y) 



yO 

'dy~^. 



hLix)-2x I'dy ^ 



(37) 
(38) 



e{x)^e{x) + --^, 



(39) 



/t(x) 



/t(x) 



y 



(40) 
(41) 



h{x) = 
ht^'\x) 



h(x) -2x i\y^ 

y 



dy 

y 



(42) 
(43) 



6^(2;) = ei,(x). 



(44) 



Note that often the combinations of tilde functions between brackets are denoted by a single 'interaction- 
dependent' function. 

In order to study the evolution of these functions, we consider the moments [/]„ = J dx x"^^ f{x)^ giving 



1 n ~ 1 m n — 1 

[9T\n = - [gi\n 4 [gTln + 77 [hl\n-l , 

n n M n 



li/lTJ" " 



([.9l]n+l - [ffT]n+l - ^ [hl]n ] , 



(45) 
(46) 



[hL]n 



[h 



Ml)] 

IL J 



2 n-l- mn-1 

n+1 n+1 

= -^^-p^ ([/«l]n+l - [ft-L]n+l - [gi]n 



(47) 
(48) 



In — r^\n 



M 



[fi 



71 



(49) 
(50) 
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[/it ]" 



n + 1 



[h] 
[ht 



-L(l)l _ 1 



n + 2 



[h]n+l , 



\eL\n = BL 



(51) 

(52) 
(53) 

(54) 



Actually, we need not consider the five T-odd functions separately. They can be simply considered as imaginary 
parts of other functions, when we allow complex functions. In particular one can expand the correlation functions 
into matrices in Dirac space [ pO[ | to show that the relevant combinations are {gir — i fir) which we can treat 
together as one complex function giT- Similarly we can use complex functions {h^j^ + i h^) —> h^j^, {gx + i fr) 
9Ti {hL +ih) ^ Hl, {e + i e^) e. 



V. EVOLUTION EQUATIONS 



In case of autonomous evolution of a function / one has 



dT 27r y 



(55) 



where r = InQ^ and P^^^ are the splitting functions. Using moments Alf of these splitting functions - the 
anomalous dimensions -, this is 



iifUr)-^Ain [/]„(.). 



(56) 



This a pplies to the leading order results for the non-singlet twist-2 functions (with the usual + prescription) 



2 ^ ' (1-/3)^ 



(57) 
(58) 



and the large Nc result in leading order for the interaction-dependent functions p3| 



p[/T(/3) = 



with c = —1 for c/T, c = —3 for and c = +1 for e. The corresponding anomalous dimensions are 

3 



(59) 



, -1 2yl 

2 n{n + 1) ^ j 



" 1 



(60) 
(61) 



and for the interaction-dependent functions in the large Nc limit 



aISt] _ ^ 



2 + n-2E- 



(62) 



7 



2 n ^7 



2 



1 3 

- H 2 > - 

2 n ^ j 



(63) 
(64) 



Using the moment analysis of the previous section, it is straightforward to find that the evolution of gjj. is 
driven not only by this function itself but also by a higher moment of gi and a similar situation for hfj^^\ In 
the large Nc limit (Cf Nc/2) one obtains (omitting mass terms) 



An 



12 n + 1 ^ j 



L^It]" + bl]n+l , 



«s(t) 

47r 



1 

2 ~ n" 



n 

--2yi 



[h 



J-(l)l 
IL J" 



n + 1 



[/il]n+l 



(65) 
(66) 



or in terms of the functions of light-cone momentum fractions 



dT 



as(r) 
47r 



N, 



dy 



1 +xy 
-6{y~x) + — — 

2 y^{y-x) + 

1 3x2 - 
-d(y-x) + — — 

2 y^{y-x)+ 



(1)/ N ^ I N 

9iT[y,v + —9i[y,T] 

V 

y 



(67) 
(68) 



Next we note that apart from a 75 matrix the operator structures of the T-odd functions /j^^'^'' and hf^^"^ are in 
fact the same as those of g[}^ and hfj^^'' (they can be considered as the imaginary part of these functions |30|). 
This implies that for the non-singlet functions, one immediately can obtain the evolution of the T-odd functions, 



^j. U IT J" ~ ■ 



An 

«s(t) 
47r 



2 n+1 p-^j 



-L(l)i 
IT J"' 



1 

2 n" 



n 



(69) 



(70) 



Furthermore, for the chiral-odd functions, which do not mix with a gluon distribution, there is no difference 
between the non-singlet and the singlet evolution. 

In the large Nc limit, the evolution equations for the non-singlct T-odd functions are of simple diagonal form 
with splitting functions 



r 

pifiT 'l(/3) = 



Nc 
2 

Nc 
2 



2'('-^)+(W)7 
2'('-^)+(W)7 



(71) 
(72) 



Actually, we also obtain the anomalous dimensions (and splitting functions) of the T-odd twist-3 functions using 



VI. FRAGMENTATION FUNCTIONS 



Just as for the distribution functions one can perform an analysis of the soft part describing quark fragmen- 
tation. One needs § 
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Ay(z,fc,) - / — — 



(27r)3 



rr(0|Z^(oo, ^'.(Ol^'/., ^)(^/., ^IV-.-WZ^lO, oo)|0) 



(73) 



Note that because of the definition of the hght-hke vectors n± (via P and Ph), the role of these vectors for 
fragmentation functions wiU be interchanged with respect to the distribution functions. For the production of 
unpolarized or spin-1/2 hadrons h in semi-inclusive hard scattering processes one needs to leading order in l/Q 
the correlation function iflOi 



A(z, kr) = zDi{z, k'^) fi_ + zD^j,{z, k 



I2\ '^A'fpo' " '^-'^T^hT 



zGis{z, k'^) ■fi_'j5 + zHit{z, k'^ ) 



Mh 



■izH^{z,k'^) 



2Mh 



O 



P.- 



(74) 



We used the shorthand notation 

Gisiz,k^) = Sh, G,L{z,k'^) 



I2\ I i^T ' ^hr) 



GiT{z,k'^), 



(75) 



etc. The arguments of the fragmentation functions are z = Pf~/k~ and k'^ = —zk^- The first is the 
(lightcone) momentum fraction of the produced hadron, the second is the transverse momentum of the pro- 
duced hadron with respect to the quark. The fcr-integrated results are, using F{z) = J (Pk!j, F{z, k'^) and 
i^(")(z) = Jd^K {kl/'iMlY F{z,k'^) = Jd^K {k'^/2z^MfXF{z,k'^), 



A*""*"'^(z) 
Ag(z) 



fi_ + ShL 75 H 75PhT ) 



z z 
Mh / E(z) , Gt{z) 

Pu 



Hl{z) I 

l5PhT + ^hL 75 - 



E^iz) Dt{z) „ ^.H{z) [fi. 

-I ShL 75 4 Ip^hTa + « 



ShT - ShL 



Dt^\z) 



z 2 

-(1) 



Mi/p7 ^+^hT * 



Ht' \z) [7",7i_ 



(76) 



(77) 



(78) 



In the twist-3 functions one can again isolate the interaction-dependent parts as done for the distribution 
functions. They are now given by 



Gt{z) = Gt{z) - zGW(z) -^zHiiz), 

Mh 

Hl{z) = Hl{z) + 2zHtP{z) - zG,{z), 

J-VJ-h 

E{z) = E{z)~^zD,{z), 
Mh 

Dt{z)^ DT{z) + zDf^^\z), 
H{z) ^ H{z) + 2zH^^^\z), 
El{z)^El{z). 



(79) 

(80) 

(81) 

(82) 
(83) 

(84) 



For the fer-integrated or the fc^/2Af/i- weighted fragmentation functions all results are obtained from the dis- 
tribution functions by replacing x — > 1/z and f...{^) — *■ F)(zji/z, g...ix) G,,,(z)/z and h,,,{x) — > H,,Xz)/z. 
The same applies to the relations from Lorcntz invariance |1GJ34[] 



Gt{z)^Gi{z)-z^ — 



G«(z) 



(85) 
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HL{z) = Hiiz)+z^ — 



Dt{z) = r 



dz 



z 



Expressing the functions in twist-2 functions and interaction-dependent functions gives 



z J 


dy 

Z 


Gi{y) 

y2 


m 


G^k^) = - J 


f dy 

Z 


Gi{y) 


m 


Hl{z) = -2^ 


t dy 

Z 


Hi{y) 
y 


m 



Hi 



J z 



dy 



Hi{y) 



Gt{z) 



I 



y2 



Hi{y) 



GT{y) 



dy 2 

y 



2Gi(z) + 2 



j\yGi{y) + Hl{z) + 2J^ 



zH. 



IL 



\z) = [ dy 

J z 



y Mh 



dy Gi {y) - I dy 



+ 

HL{y) 
y 



H,{z) + 2l\y^ 

y 



(86) 
(87) 
(88) 

(89) 
(90) 

(91) 
(92) 



E{z)=E{z) + —zDi{z), 



(93) 



Driz) 



D-^\z) 



Dt{z) 



J z 



y2 



(94) 

(95) 



H{z) 
zHt^%) = -f 

J z 

El{z) = El{z). 



H{z) +2 l\y ^^^^ 



y 



dy 



H{y) 



(96) 
(97) 

(98) 



The relations for the moments of fragmentation functions can be obtained from the above equations or from 
the results of the distribution functions via the replacements n — > —n followed by [f]-n — * [D/z]n = [D]n-i. 
This yields 



[Grin = — [Gl]n H ^ [GT\n + JJ ^ [Hi]n+1 

n + 1 n + 1 Mh n + 1 



Mh n- 

^ ([Gi]„ - [GtU - ^ [HiUi I , 



(99) 
(100) 



[HL]n = [Hi]n H ^— [HL]n + TT — [Gl]„+1 

n n Mh n 



Wth^+l = - ( [Hl]n - [HL\n - [Gi]„+i I , 



(101) 
(102) 
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(103) 



n + 1 



1 



71+1 



[DtU 



(104) 
(105) 



[H]n = ^ [H]n , 



(106) 
(107) 



[£^L]n = [EL]n ■ 

The autonomous evolution equations are again of the form 

or via the (usual) moments — dz p[^l(z) of the splitting functions, 



^4^1 [DUr) 

ZTT 



(108) 



(109) 



(110) 



For the leading order contributions the analytic structure of the corrections for fragmentation functions is similar 
as for distribution functions. We note a (generalized) Gribov-Lipatov reciprocity, summarized by the following 
procedure. The splitting functions for distribution functions /(x,t) and corresponding fragmentation functions 
z D{z, t) are related by 



p[-^l(/3) 



(1-/3)+' 



(Ill) 

(112) 



This relation works for the twist-2 fragmentation functions and the interaction-dependent functions ||2^, for 
7V(/3) being (at most a quadratic) polynomial in f3. In the case of the twist-2 functions the functional form of 
the splitting functions is the same for distribution and fragmentation functions. This is no longer true for the 
interaction-dependent functions. For the anomalous dimensions of distribution and fragmentation functions the 
relation becomes 



1 

AlP = A{n) - 2y - = A{n) - 2-/e - 2^{n + 1), 
1 

n 

43i = A {-{n + 1)) - 2^E - mn +l)=A {-{n + 1)) - 2 ^ - 



(113) 
(114) 



where A{n) is a rational function. We have not yet investigated the wider applicability of the above relations. 
We find for the twist-2 fragmentation functions the familiar results, which obey the original Gribov-Lipatov 
reciprocity relation aK^ — A\^J^ between the twist-2 distribution functions / = fi,gi,hi and fragmentation 
functions D = Di,Gi, Hi, 



aIDi] _ a[Gi] _ ^ 
— ^n+1 — ^'-F 



1 



n{n + 1) 



1 



A, 



Hi] 
n+l 



2 ^ j 



(115) 
(116) 
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In the large Nc limit, our generalized reciprocity relations in Eqs. (113) and (114) applied to Eqs. ( |62| ) - ( |64| ) 
give the results for the interaction-dependent functions [^sj^ 



a[Gt] 



aIHl] 



2 
2 



2 



1 



n + 1 
1 

n+1 
3 

n+1 



1 

n 
n 



(117) 
(118) 
(119) 



Again one then also knows ^[^^1 = ^[^^1 = A^^^\ A^"^ = A^^^ = A^"^^ and A^'^^^ = A^^^^ = Al^l = A^^\ 

Using the moment analysis (the reciprocity relations cannot be used straightforwardly) one obtains, omitting 
the mass terms. 



An 



d ±(1) _ a.(T) 




{n - l){n + 1) 



[Gl]n } , 



(120) 
(121) 



with in this case mixing with a lower moment of the twist-2 functions. In terms of the functions of lightcone 
momentum fractions one finds 



dy 

1 



dy 



t: o(y ~ z) + — 

1 X , iy~ z 

-d(y - z) + — — 

2 y{y~z)+ 



yG[%,r)-y^^G,{y,r)}, 

yH^L \y^ V + 7. ^1 (y' V 



y 



(122) 
(123) 



Given the fact that, apart from an additional 75, the operator structure for the T-odd Sivers and Collins 
functions, of^^^ and H^'"^\ are the same as those of G^^l^ and H^jp but without mixing with Gi or Hi, one 
finds in the large limit an autonomous evolution for the T-odd functions, with anomalous dimensions 

1 



A 



n+1 



A. 



n+1 



2 


1 

2^ 


1 

n 


n 










Nc 
2 


1 

2 ~ 


3 

h - 
n 


n 



(124) 
(125) 



corresponding to splitting functions 

r r-.^(l)l 
p[zD,^ l(/3) 



2 

Nc 
2 



(126) 
(127) 



The results Eqs. (125) and (127) should prove useful for studies of the Collins effect and Eqs. (124) and ( |126| ) 
for studies of transversely polarized A production [p5[. 
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TABLE I. The anomalous dimensions from which the large behavior of the moments, proportional to [Qji(Q^)]''", 
is obtained. Defining the moments a„ taking out the factor Cf or from the anomalous dimensions yl„, one has for 

the twist-2 functions d„ = — 2a„ Cf /bo with bo = (11 A'^c — 2 Nf)/3, while for the large Nc results one has d„ = — 3 a„/ll. 
Also indicated is the charge conjugation behavior of the functions, f{x) — ±/(— s). 



function 


c 




CL2 


a3 


validity 


J -L 







-7/6 


-25/12 











-7/6 


-25/12 




hi 




-1/2 


-3/2 


-13/6 




qt and fT 




-1/2 


-2 


-17/6 


larffe JV^ 


Hl and h 




-5/2 


-3 


-7/2 


o ' c 


e 


+ 


-1-3/2 


-1 


-13/6 


large Ac 


q\V, and f."^^^^ 




-2 


-17/6 


-41/12 


large A^c 


and ftf 


+ 


-3 


-7/2 


-47/12 


large Nc 


zDi 







-7/6 


-25/12 




zGi 


+ 





-7/6 


-25/12 




zHi 




-1/2 


-3/2 


-13/6 




zGt and zDt 


+ 


-2 


-17/6 


-41/12 


large Ac 


zHl and z_ff 




-1 


-13/6 


-35/12 


large Ac 




+ 


-3 


-7/2 


-47/12 


large A^c 


zG^Y^ and zOf^^ 




-1/2 


-2 


-17/6 


large A^c 


zHt-^'^ and 


+ 


+3/2 


-1 


-13/6 


large Ac 



VII. DISCUSSION AND CONCLUSIONS 

Our goal was to obtain the evolution equations of the functions that appear in azimuthal spin asymmetries. 
These pr-dependent functions appear in asymmetries that are not suppressed by explicit powers of the hard 
momentum. But as functions of transverse momentum they are not of definite twist, which implies that in 
order to obtain the evolution equations one has to calculate corrections to higher twist operators as well. 
For the first p|/2M^ moment (transverse moment) of these pr-dependent functions, such as for the Collins 
fragmentation function H^^^'' = J (Pk!j,k'^/2z'^M'^ H^{z,k'^), we obtain DGLAP-like evolution equations. 
Such moments appear in cross sections weighted with the momentum g", where only the directional (azimuthal) 
dependence remains. For explicit examples we refer to Refs. ||l^,|ll|. In case one does not weight the transverse 
momentum integration of the differential cross section, one is only sensitive to the twist-2 functions fi,gi and 
hi (and their fragmentation counterparts), but in case one weighs with one or more powers of the observed 
transverse momentum, one becomes sensitive to the functions hf^"^\ fiT^\ h^'"^^ (and their fragmentation 
counterparts), which are functions of the lightcone momentum fraction x only. 

In the large- A^c limit, the non-singlet evolution of these functions involves only the functions themselves and 
(in the T-even case) only well-known twist-2 functions. For the chiral-odd functions the equations also apply to 
the singlet case, since there is no mixing with gluon distribution functions. The large- A'^j, evolution equations 
are expected to be good approximations to the full evolution equations which are not of this simple form (cf. 
Ref. [Q), because of the appearance of two-argument twist-3 functions as in Eq. (p^). It is not excluded that 
the first 1/Nc correction to the result obtained here may still lead to autonomous evolution equations, but we 
will not address this issue here. Especially the (large Nc) evolution equation we have obtained for H^'^^\ 

yHt^'\y,r), (128) 

should prove useful for the comparison of data on Collins function asymmetries from different experiments, 
performed at different energies. 

It is worth investigating the large Q behavior of the solutions to the various evolution equations. For this 
purpose we have given the first 3 anomalous dimensions for the different functions in Table ||. First we note that 
all (diagonal) anomalous dimensions oi g[l^ , hf^^\ ff^^^ and hf^-^'^ are negative, implying that these functions 
will vanish asymptotically {Q^ —> oo), except that for the T-even functions there is mixing with gi and /ii, but 
this does not alter the conclusion. 



Ml) 



47r 



Nr 



dy 



-6{y~z) + — — 

2 y(y~z)+_ 



13 



For the fragmentation counterparts the conclusion is similar, except for the fact that the lowest anomalous 
dimensions of zH^^^'^ and zH^^^"^ are positive, potentially leading to divergent behavior of the functions as 
— > oo. However, here we recall the Schafer-Teryaev sum rule p% 

J dz zHI^^\z) = 0, (129) 

which states that the first moment of zH^^^"^ is zero, making the sign of the first anomalous dimension irrelevant. 
Similar sum rules hold for the other first transverse moments of fragmentation functions js^js^. All higher 
moments will vanish asymptotically. The behavior of the sum rule for the first moment of the function e is 
discussed in Ref . (D . 

In conclusion, using Lorentz invariance and the QCD equations of motion, the operator structure of the 
transverse moments of pT-dependent quark distribution and fragmentation functions can be found in terms of 
twist-2 and twist-3 operators. Knowing their, for large Nc simple, evolution one also knows the evolution of 
azimuthal asymmetries in semi-inclusive hard scattering processes. 
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